Band structure from random interactions 
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The anharmonic vibrator and rotor regions in nuclei are investigated in the framework of the 
interacting boson model using an ensemble of random one- and two-body interactions. We find 
a predominance of = 0^ ground states, as well as strong evidence for the occurrence of both 
vibrational and rotational band structures. This remarkable result suggests that such band struc- 
tures represent a far more general (robust) property of the collective model space than is generally 
0> ■ thought. 
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^ ' A recent analysis of experimental energy systematics of medium and heavy even-even nuclei suggests a tripartite 
O > classification of nuclear structure into seniority, anharmonic vibrator and rotor regions [Qj^. Plots of the excitation 
energies of the yrast states with = 4+ against — 2+ show a characteristic slope for each region: 1.00, 2.00 
and 3.33, respectively. In each of these three regimes, the energy systematics is extremely robust. Moreover, the 
transitions between different regions occur very rapidly, typically with the addition or removal of only one or two 
pairs of nucleons. The transition between the seniority region (either semimagic or nearly semimagic nuclei) and 
the anharmonic vibrator regime (either vibrational or 7 soft nuclei) was addressed in a simple schematic shell model 
^ , calculation and attributed to the proton-neutron interaction ||^ . The empirical characteristics of the collective regime 
which consists of the anharmonic vibrator and the rotor regions, as well as the transition between them, have been 
studied [QH] in the framework of the interacting boson model (IBM) Q . An analysis of phase transitions in the IBM 
has shown that the collective region is characterized by two basic phases (spherical and deformed) with a sharp 
transition region, rather than a gradual softening which is traditionally associated with the onset of deformation in 
ON ■ nuclei |@. 

- In a separate development, the characteristics of low-energy spectra of many-body even-even nuclear systems have 



i -^H ' been studied recently in the context of the nuclear shell model with random two-body interactions 1 11 1^ . Despite the 
I , random nature of the interactions, the low-lying spectra still show surprisingly regular features, such as a predominance 
^ of = 0"*" ground states separated by an energy gap from the excited states, and the evidence of phonon vibrations. 
^ The occurrence of these pairing effects cannot be explained by the time-reversal symmetry of the random interactions 
. A subsequent analysis of the pair transfer amplitudes has shown that pairing is a robust feature of the general 
J> , two-body nature of shell model interactions and the structure of the model space [Q. On the other hand, no evidence 
\^ ' was found for rotational band structures. 

, The existence of robust features in the low-lying spectra of medium and heavy even-even nuclei ||l|,|^ suggests that 
' there exists an underlying simplicity of low-energy nuclear structure never before appreciated. In order to address 
this point we carry out a study of the systematics of collective levels in the framework of the IBM with random 
interactions. In an analysis of energies and quadrupole transitions we show that despite the random nature (both in 
size and sign) of the interaction terms, regular features characteristic of the anharmonic vibrator and rotor regions 
emerge. Our results imply that these features are, to a certain extent, independent of the specific character of the 
interaction, and probably arise from the two-body nature of the Hamiltonian and the structure of the collective model 
space. 

In the IBM, collective nuclei are described as a system of N interacting monopole and quadrupole bosons. We 
consider the most general one- and two-body IBM Hamiltonian H — Hi + H2- The one- and two-body matrix 
elements are chosen independently using a Gaussian distribution of random numbers with zero mean and variances 



(Hip,,) ^—^v'{l + 6^p,). (1) 



Since the matrix elements of Hi and H2 are proportional to N and N{N — 1), respectively, we have introduced a 
relative scaling between the one- and two-body interaction terms of i/{N — 1). The coefficient is independent of 
the angular momentum and represents an overall energy scale. The ensemble defined by Eq. (pi) is similar, but not 
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identical, to the Two-Body Random Ensemble of iQ. In all calculations we take A'' = 16 bosons and 1000 runs. For 
each set of randomly generated one- and two-body matrix elements we calculate the entire energy spectrum and the 
B{E2) values between the yrast states. 

Just as in the case of the nuclear shell model we find a predominance (63.4 %) of = 0"'" ground states; in 
13.8 % of the cases the ground state has — 2+, and in 16.7% it has the maximum value of the angular momentum 
= 32+. For the cases with a = 0+ ground state we have calculated the probability distribution of the energy 
ratio R = [E{4+) - E{0+)]/[E{2+) - E{0+)]. Fig. |l| shows a remarkable result: the probability distribution P{R) 
has two very pronounced peaks, one at i? ~ 1.95 and a narrower one at i? ~ 3.35. These values correspond almost 
exactly to the harmonic vibrator and rotor values (see the results for the C/(5) and SU(S) limits in Table |). No such 
peak is observed for the 7 unstable or deformed oscillator case (5*0(6) limit). 

Energies by themselves are not sufficient to decide whether or not there exists band structure. Levels belonging to 
a collective band are connected by strong electromagnetic transitions. In Fig. ^ we show a correlation plot between 
the ratio of B{E2) values for the 4+ — > 2+ and 2+ — > 0+ transitions and the energy ratio R. For the B{E2) values we 
use the quadrupole operator 

Q^(x) = (.stJ+dt42)+^(rftJ)(2) ^ (2) 

with X = ^\/7/2. For completeness, in Table | we show the results for the three symmetry limits of the IBM [||. 
In the large N limit, the ratio of B{E2) values is 2 for the harmonic oscillator (J7(5) limit) and 10/7 both for the 
deformed oscillator {SO{Q) limit) and the rotor {SU{i) limit). There is a strong correlation between the first peak 
in the energy ratio and the vibrator value for the ratio of B{E2) values (the concentration of points in this region 
corresponds to about 50 % of all cases), as well as for the second peak and the rotor value (about 25 % of all cases). 
For the region 2.3- i?- 3.0 one can see a concentration of points around the value 1.4, which reflects the transition 
between the deformed oscillator and the rotor limits (see Table |) . Calculations for different values of the number of 
bosons N show the same results. 

Despite the randomness of the interactions these results constitute strong evidence for the occurrence of both 
vibrational and rotational band structures. We have repeated the calculations for different values of the number 
of bosons N and find the same results. Since the results presented in Figs. |l| and || were obtained with random 
interactions, with no restriction on the sign or size of the one- and two-body matrix elements, it is of interest to 
compare them with a calculation in which the parameters are restricted to the 'physically' allowed region. To this 
end we consider the consistent Q-formulation [ p^ which uses the same form for the quadrupole operator, Eq. (^, i.e. 
with the same value of x, for the E2 operator and the Hamiltonian 

H = ehd-KQix)-Q{x) ■ (3) 

The parameters e and k are restricted to be positive, whereas x can be either positive or negative —V7 /2 < x < V7/2. 
The properties of the Hamiltonian of Eq. (|^) can be investigated by taking the scaled parameters rj = e/[e + 4:K{N— 1)] 
and X = 2x/'\/7 randomly on the intervals < 77 < 1 and — 1 < x ^ 1 (these coefficients have been used as control 
parameters in a study of phase transitions in the IBM H,^). In Figs. ^ and ^ we show the corresponding probability 
distribution and correlation plot for the consistent Q-formulation of the IBM with realistic interactions. Although 
in this case the points are concentrated in a smaller region of the plot than before, the results show the same 
quantitative behavior as for the IBM with random one- and two-body interactions. In Fig. ^ we have identified each 
of the dynamical symmetries of the IBM (and the transitions between them). There is a large overlap between the 
regions with the highest concentration of points in Figs. || and ^. 

In conclusion, we have studied the IBM using random ensembles of one- and two-body Hamiltonians. It was found 
that despite the randomness of the interactions the ground state has = 0+ in 63.4 % of the cases. For this subset, 
the analysis of both energies and quadrupole transitions shows strong evidence for the occurrence of both vibrational 
and rotational band structure. These features arise from a much wider class of Hamiltonians than are generally 
considered to be 'realistic'. This suggests that these band structures arise, at least in part, as a consequence of the 
one- and two-body nature of the interactions and the structure of the collective model space, and hence represent 
a far more general and robust property of collective Hamiltonians than is commonly thought. This is in qualitative 
agreement with the empirical observations of robust features in the low-lying spectra of medium and heavy even-even 
nuclei 

A similar situation has been observed in the context of the nuclear shell model with respect to the pairing properties 
0,0 which were formerly exclusively attributed to the particular form of the nucleon-nucleon force. On the other 
hand, the random IBM Hamiltonians studied in this Letter display not only vibrational-like phonon collectivity but, in 
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contrast to the results in |T^,|l^, also imply the emergence of rotational bands. The IBM is based on the assumption 
that low-lying collective excitations in nuclei can be described as a system of interacting monopole and quadrupole 
bosons, which in turn are associated with generalized pairs of like-nucleons with angular momentum L = and 
L = 2. It would be very interesting to establish whether rotational features can also arise from ensembles of random 
interactions in the nuclear shell model, if appropriate (minimal) restrictions are imposed on the parameter space. 

It is a pleasure to thank Stuart Pittel and Rick Casten for interesting discussions. This work was supported in part 
by DGAPA-UNAM under project IN101997. 
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TABLE I. Energies of B{E2) values in the dynamical symmetry limits of the IBM In the f/(5) and SO(6) limits we 
show the result for the leading order contribution to the rotational spectra. 
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FIG. 1. Probability distribution P{R) of the energy ratio R = [E{4+) - E(0+)]/[E{2+) - E{0+)] with / P(R)dR = 1 in the 
IBM with random one- and two-body interactions. 
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FIG. 2. Correlation between ratios of B{E2) values and energies in the IBM with random one- and two-body interactions. 
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FIG. 3. As Fig. |l[ but in the consistent Q-formulation of the IBM. 
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FIG. 4. As Fig. pi but in the consistent Q-formulation of the IBM. 



